Graphene quantum dots embedded in hexagonal boron nitride sheets 
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We have carried out first-principles calculations on electronic properties of graphene 
quantum dots embedded in hexagonal boron nitride monolayer sheets. The calcula- 
tions with density functional theory show that the band gaps of quantum dots are 
determined by the quantum confinement effects and the hybridization of ir orbitals 
from B, N and C atoms. The energy states near the Fermi level are found to be 
strongly localized within and in the vicinity of the quantum dots. 
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Graphene has attracted a great deal of interest since it was exfoliated from bulk 
graphite.- 1,2 Because graphene has a number of unusual properties such as massless fermions 
and high carrier mobility, graphene materials are expected to play an important role in fun- 
damental research and in future electronic device applications. 3 ' 4 Two-dimensional graphene 
has zero band gap with linear energy dispersion near the Fermi level. Many efforts have 
been devoted to opening a band gap in order to fit various needs. One promising method is 
to confine charge carriers in two dimensions by cutting graphene into ribbons," of which 
the zigzag graphene nanoribbons are of special interest because of the spin-polarized edge 
states predicted to be half-metallic upon applying transverse electric field.- 

Electronic states in graphene could be further confined in the basal plane to make quan- 
tum dots (QDs) and there are already a few investigations on this idea. Zhang et al. 10 
showed that the surface corrugation of graphene on Ru(0001) can form QDs. Ponomarenko 
et al— carved out graphene dot devices behaving as single-electron transistors and exhibit- 
ing coulomb blockade behavior. Yan et al— found that solution-processable graphene dot 
can be used as light absorbers for photovoltaic applications. The QDs also could be hosted 
in other two-dimensional systems. Singh et al— using density functional theory and the 
tight-binding method, proposed that H-vacancy clusters in graphane could form graphene 
quantum dots having size dependence of energy gaps typical for confined Dirac fermions. 

Along with the development of graphene, rapid progress has been made on the study of 
boron nitride materials. Recently graphene-like hexagonal boron nitride (h-BN) sheets have 
been fabricated experimentally.— - — When C atoms mix with boron nitride, many interesting 
properties can emerge. Even without the need of externally applying electric field, hybrid C- 
BN nanotubes and nanoribbons could exhibit intrinsic half-metallicity— ~— and the C doping 
can induce magnetism in BN nanotubes.— 

In the present letter, we report the first-principles calculations on the graphene QDs with 
h-BN sheets as host materials. Our results show that the hybridization between 7r orbitals 
of B, N and C atoms determines the energy gap together with the quantum confinement 
effect. The states close to Fermi level Ep are found strongly localized around the C-doped 
region. 

The calculations were performed with spin-unrestricted density functional theory as im- 
plemented in the SIESTA code.— The norm-conserving Troullier- Mart ins pseudopotentials 
were used to model the interaction between ionic cores and valence electrons and double- 
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C basis set plus polarization orbitals were employed to construct the valence orbitals. The 
exchange- correlation functional was described with Per dew- Burke- Ernzerhof approximation. 
A grid cut-off of 200 Ry was used. The graphene QDs embedded in h-BN monolayer are 
simulated by doping C atoms in a h-BN sheet of rectangular shape. A vacuum region of 
16.5 A was added along the direction normal to the sheet plane to avoid the interaction 
between periodic images of supercells. Relaxations were carried out on all structures inves- 
tigated until the force on each atom is less than 0.04 eV/A. The k-points sampling schemes 
of lxlxl and 2x2x1 were employed for geometry optimization and density of states (DOS) 
analysis, respectively. 

The QDs of high symmetry, as depicted in Fig. [TJ were chosen in this study. They were 
constructed by replacing all B and N atoms enclosed within a circle of diameter d by C 
atoms with the total number of atoms fixed at 512. We varied d from 2.85 A to 12.55 A, 
corresponding to the total C atoms n = 6, 12, 24, 36 and 54, to model QDs of different sizes 
denoted by C n . 



FIG. 1. Structural model of QD C36 embedded in h-BN sheet. B, N and C atoms are represented 
by gray, blue and dark-green spheres, respectively. 

The stability of QDs is estimated from the formation energy which can be expressed as 

E f = E tot - m/x BN + 2m/i C , 

where i?tot is the total energy of QDs and anc ^ /^C are ^ ne chemical potentials of BN 
pair and C atom, respectively. /igN an< ^ A*C are calculated from the cohesive energies of 
a BN pair in an infinite h-BN monolayer and of each C atom in a single graphene sheet, 
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respectively. Here m represents the total number of BN pairs replaced by C atoms. We 
display the formation energies per C atom as a function of QDs size in Fig. [2j Roughly, the 
formation energies decrease with increasing quantum dot size. The detailed size dependence, 
however, is oscillating instead of monotonic. Note that the QDs Cg, C24, and C54 comprised 
of complete aromatic rings have lower formation energy than C12 and C36. This shape 
preference was also observed in graphene nanoribbons embedded in boron nitride sheets^ 
and graphene QDs originating from H- vacancies in graphane. 13 




FIG. 2. Formation energy Ef as a function of QDs diameter d. 

For pure h-BN sheet, our calculated energy gap is 4.57 eV comparable to 4.61 eV reported 
from calculations with plane wave basis.— After pairs of B and N were substituted by C, 
the energy gap is reduced. Fig. [3] depicts the change in energy gap for various values of d. 
When increasing the QDs size, the energy gap reduces from 3.57 eV to 1.64 eV. In isolated 
graphene flakes terminated with H atoms, the gap change with size is mainly determined 
by quantum confinement effect since H atom only has one s electron forming sp 2 a bonds 
far away from the Fermi level. This is similar to the QDs from graphane except the sp 3 a 
bonds are formed instead. When the QDs are surrounded by B and N atoms in the present 
study, the situation is different. In addition to the confinement effect, the 2p electrons of 
B and N will form n bonding near the Fermi level with 2p electrons of C. We can see this 
point based on the DOS analysis. 

In Fig. H] we show the total density of states (TDOS) and partial density of states (PDOS) 
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FIG. 3. Energy gap E g as a function of QD diameter d. Solid cure is drawn as a guide to the eyes. 



decomposed onto 2p orbitals for h-BN sheet and QD C36 {d = 10.50 A). For pure h-BN sheet, 
the top of valence bands are dominated by 2p orbitals from N and the bottom of conduction 
bands mainly originate from 2p orbitals of B atoms. The introduction of graphene QDs has 
almost no effect in the energy range of valence and conduction bands of h-BN sheet. However, 
the near-Ep DOS are largely modified and several peaks show up. The peaks 0.85 eV below 
and above Ep are almost dispersionless and there is also strong orbital hybridization between 
2p orbitals of B, N and C atoms. The mixing of 2p orbitals will form bonding and anti- 
bonding states. The energies of the 2p orbital are different with B higher than N and C 
in the middle. So near the Fermi level, only 7r c . B and 

^C-N 

the relevant bands. To see 

how these states distribute in real space, we show the local density of states (LDOS) for 



the highest occupied molecular orbitals (HOMO) in Fig. 5(a) and for the lowest unoccupied 



molecular orbitals (LUMO) in Fig. 5(b) These energy states are strongly localized in the C 
region and also have considerable distribution on B and N atoms at the interface of C-BN 
in the form of 7r c _ B and n* _ N states. In the interior region, compared to HOMO, the LUMO 
has more nodes which increase the kinetic contribution in energy. The discussion based on 
C36 also applies on other QDs in which we observed similar orbital hybridization and carrier 
localization. 
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FIG. 4. (a) TDOS and (b) PDOS of h-BN sheet, (c) TDOS and (d) PDOS of QD C 36 - PDOS 
represent 2p orbital contributions of B, N and C as indicated. The vertical dotted line indicates 
the Fermi level as zero energy reference. We apply a broadening of 50 meV to the energy levels. 
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In summary, we have performed pseudopential calculations with numerical atomic orbitals 
basis to investigate the electronic properties of graphene QDs embedded in h-BN sheet. 
Our results show that the orbital hybridization between 2p orbitals of B, N and C and 
the quantum confinement together determine the energy gaps of QDs. The dispersionless 
energy states near Fermi level are strongly confined in the QDs region. Great progress in 
the synthesis of single layer h-BN sheets and graphene-BN composites 27 makes the kind 
of quantum dots discussed in our work potentially promising candidates in electronic and 
optoelectronic applications. 

We gratefully acknowledge research support from the NSF and NRI through the Brown 
University MRSEC program and the NSF through the grants CMMI-0825771, DMS-0914648 
and DMS-0854919. 



8 



REFERENCES 



1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. 



Grigorieva, and A. A. Firsov, |Science 306, 66 6 (2004). 

2 K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. Khotkevich, S. V. Morozov, and 
A. K. Geim, Proc. Natl. Acad. Sci. 102, 10451 (2005). 

3 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, 
S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005). 

4 Y.-M. Lin, C. Dimitrakopoulos, K. A. Jenkins, D. B. Farmer, H.-Y. Chiu, A. Grill, and 
P. Avouris, |Science 327, 662 (2010)[ 

5 M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, |Phys. Rev. Lett. 98, 206805 (2007)| 
6 J. Campos-Delgado, J. M. Romo-Herrera, X. Jia, D. A. Cullen, H. Muramatsu, Y. A. Kim, 
T. Hayashi, Z. Ren, D. J. Smith, Y. Okuno, T. Ohba, H. Kanoh, K. Kaneko, M. Endo, 



H. Terrones, M. S. Dresselhaus, and M. Terrones, Nano Lett 8, 2773 (2008) 



7 X. Jia, M. Hofmann, V. Meunier, B. G. Sumpter, J. Campos-Delgado, J. M. Romo- 
Herrera, H. Son, Y.-P. Hsieh, A. Reina, J. Kong, M. Terrones, and M. S. Dresselhaus, 
Science 323, 1701 (2009) . 



8 L. Jiao, L. Zhang, X. Wang, G. Diankov, and H. Dai, Nature 458, 877 (2009). 
9 Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature 444, 347 (2006). 

10 H. G. Zhang, H. Hu, Y. Pan, J. H. Mao, M. Gao, H. M. Guo, S. X. Du, T. Greber, and 



H.-J. Gao, |J. Phys.: Condens. Matter 22, 302001 (20 10), 
n L. A. Ponomarenko, F. Schedin, M. I. Katsnelson, R. Yang, E. W. Hill, K. S. Novoselov, 
and A. K. Geim, |Science 320, 356 (2008) 



12 X. Yan, X. Cui, B. Li, and L.-S. Li, |Nano Lett. 10, 1869 (2010)] 

13 A. K. Singh, E. S. Penev, and B. I. Yakobson, ACS Nano 4, 3510 (2010) 



14 D. Pacile, J. C. Meyer, Q. O. Girit, and A. Zettl, |Appl. Phys. Lett. 92, 133107 (2 008). 
15 W.-Q. Han, L. Wu, Y. Zhu, K. Watanabe, and T. Taniguchi, 

|Appl. Phys. Lett. 93, 223103 (2008)| 
16 C. Jin, F. Lin, K. Suenaga, and S. Iijima, |Phys. Rev. Lett. 102, 195505 (2009)[ 
17 N. Alem, R. Erni, C. Kisielowski, M. D. Rossell, W. Gannett, and A. Zettl, 



Phys. Rev. B 80, 155425 (2009) 



is 



L. Song, L. Ci, H. Lu, P. B. Sorokin, C. Jin, J. Ni, A. G. Kvashnin, D. G. Kvashnin, 



9 



J. Lou, B. I. Yakobson, and P. M. Ajayan, |Nano Lett. 10, 3209 (2010) 
19 F. Miiller, S. Hlifner, H. Sachdev, S 



Gsell, 



and M. Schreck, 



Phys. Rev. B 82, 075405 (2010) 



20 



S. Dutta, A. K. Manna, and S. K. Pati, |Phys. Rev. Lett. 102 , 096601 (2009) 



21 



B. Huang, C. Si, H. Lee, L. Zhao, J. Wu, B.-L. Gu, and W. Duan, 



Appl. Phys. Lett. 97, 043115 (2010) 



22 J. M. Pruneda, |Phys. Rev. B 81, 161409 (2010) 
23 Y. Ding, Y. Wang, and J. Ni, |Appl. Phys. Lett. 95, 123l"05~(2009) 



24 R. Q. Wu, L. Liu, G. W. Peng, and Y. P. Feng, |Appl. Phys. Lett. 86, 122510 (2005)| 
25 J. M. Soler, E. Artacho, J. D. Gale, A. Garca, J. Junquera, P. Ordejn, and D. Snchez- 

Portal, J. Phys.: Condens. Matter 14, 2745 (2002). 
26 H. §ahin, S. Cahangirov, M. Topsakal, E. Bekaroglu, E. Akturk, R. T. Senger, and 

S. Ciraci, |Phys. Rev. B 80, 155453 (2009)| 
27 L. Ci, L. Song, C. Jin, D. Jariwala, D. Wu, Y. Li, A. Srivastava, Z. F. Wang, K. Storr, 



L. Balicas, F. Liu, and P. M. Ajayan, Nature Mater. 9, 430 (2010) 



10 



